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INTRODUCTION 


The  model  presented  in  this  technical  memorandum  was  developed  in 
support  of  the  expendable  jammer  study  currently  being  conducted  at  the 
Naval  Air  Development  Center. 

The  purpose  of  this  model  is  to  evaluate  the  relative  effectiveness 
of  friendly  jamming  of  enemy  acquisition  radars  in  terms  of  the  tracking 
errors  introduced  into  the  enemy  acquisition  radar  prior  to  target  hand 
over  to  the  enemy  fire  control  radar  operators.  The  effectiveness  of 
friendly  jamming  against  enemy  acquisition  radars  was  evaluated  and  is 
reported  on  in  reference  (a) . 

Typical  results  from  the  Radar  Entrapment  Model  developed  herein 
are  described  in  reference  (a) .  The  radar  entrapment  model  is  presently 
programmed  on  the  NAVAIRDEVCEN  (Naval  Air  Development  Center)  CDC-3200 
computer,  and  is  useful  for  the  evaluation  of  jamming  performance  in 
terms  of  the  errors  introduced  into  the  enemy  acquisition  radar  in  terms 
of  penetrator  track  prediction. 

At  present,  the  model  is  being  expanded  to  include  the  effects  of 
expendable  jammers,  that  is,  multiple  jamming  strobes,  expendable  jammer 
ballistic  trajectory  effects,  and  the  effects  of  possible  enemy  radar 
burnthrough .  • 

The  basic  radar  entrapment  model  presented  herein  assumes  that  the 
enemy  acquisition  is  tracking  a  target  in  error  and  is  estimating  the 
actual  penetration.  Track  information  is  handed  over  to  the  enemy  fire 
control  radar  which  is  required  to  position  itself  (in  space)  so  as  to 
effect  an  intercept  with  the  SAM  (Surface-to-Air  Missile)  envelope. 
However,  due  to  track  errors  introduced  by  the  jamming,  the  radar  position 
vector  is  most  likely  in  error.  The  model  evaluates  this  error  in  terms 
of  the  probability  that  the  aircraft  will  be  physically  present  in  the 
fire  control  radar  entrapment  region  when  the  fire  control  radar  is 
activated. 
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DISCUSSION 


Description 

We  assume  a  friendly  airplane  (or  a  cluster  of  airplanes)  carrying 
on  board  jamming  devices  -  describes  three  dimensional  straight  line 
constant  velocity  motion. 

Two  enemy  radar  sources  are  present.  One  source,  the  acquisition 
radar,  tracks  the  airplane  until  the  enemy  estimates,  based  on 
observed  data^  of  when  and  where  the  airplane  will  most  likely  be  first 
within  the  lethal  detection  'region  about  the  other  source,  the  fire 
control  radar.  In  order  to  seize  the  earliest  opportunity  to  fire  at 
the  airplane,  the  latter,  close  to  which  an  enemy  missile  launcher  is 
located,  then  pivots  angle -wise  about  its  fixed  location  so  that  one^ 
comer  of  its  lethal  detection  region  is  placed  at  the  estimated  posi¬ 
tion  of  entrance  of  the  airplane  within  the  effective  missile  range. 

The  lethal  detection  region  for  the  fire  control  radar  is  defined 
as  that  region  bounded  by  the  maxinial  and  minimal  effective  missile 
ranges  and  also  within  the  finite  detection  cone  of  the  radar,  the 
position  of  which  is  angular  dependent  only.  (See  figure  1.) 

Estimates  are  also  made  by  the  enemy  of  the  estimated  exit  time  of  the 
airplane,  and  time  gates  are  established  about  the  estimated  entrance 
and  exit  times  from  the  lethal  region.  These  calculations  determine 
the  turn-on  and  turn-off  times  of  the  fire  control  radar.  Since  the 
longer  the  fire  control  radar  is  turned  on,  the  more  likely  the  enemy 
will  itself  be  detected  and  fired  upon  by  the  aircraft,  an  upper  bound 
is  also  imposed  on  the  elapsed  time  between  turn-on  and  turn-off  times. 

Standard  regression  analysis  is  used  in  obtaining  the  enemy's 
estimated  path  of  the  target.  Outcomes  are  then  generated  at  various 
sample  times  and  a  determination  is  made  whether  the  airplane  will 
actually  be  in  the  lethal  region  between  enemy  turn-on  and  turn-off 
times,  for  a  given  run  (set  of  sampling  outcomes).  The  estimate  of 
probability  of  (lethal)  detection  is  then  just  the  percentage  of  runs 
in  which  the  airplane  passes  through  the  lethal  region  between  turn-on 
and  turn-off  times. 
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Model  Development 

Parameters  of  equations  of  motion  of  airplane  T: 
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Through  jamming  noise,  observations 


where  for  i  =  l,...m  (observations  are  made  by  another  radar  located 
at  another  source) : 

f.  is  the  observed  x-y  range  of  T  by  the  acquisition  radar 
at  time  ti 

O 

0.  is  the  observed  azimuth  angle  of  T  by  the  acquisition  radar 
at  time  ti 

fJ.  is  the  observed  elevation  angle  of  T  by  the  acquisition 
radar  at  time  ti 

I 

(  See  appendix  A  for  further  evaluations. ^ 

V  ,  V  ,  V  ,  X  ,  y  ,  z  are  totally  unknown  to  radar 
X  y  z  o  "^0  o 

The  relations  between  x-y-z  rectangular  coordinates  and  r-0-f} 
coordinates  are 


X  =  rcosO 
y  =  r  Sin  6 
z  =  r  t  ^  n  0 
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It  follows  from  standard  results  (reference  (a))  that  best 
estimate  of  X 

=  minimal  covariance  matrix  absolute  unbiased  estimate 

=  generalized  least  squares  estimate,  etc. 

based  on  data  x^,  y^, 
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The  above  results  constitute  the  ’'direct"  regression  solution. 
However,  a  6  by  6  matrix  must  be  inverted  for  each  new  observation 
time,  Recursive  filtering  as  given  by  the  Kalman-Buoy  filter 

reduces  this  problem; 
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Then  if  we  let  X[p+lJ  be  the  "best"  predicted  value  of  X  .  based  on 
data  up  to  time^tp,  we  can  show  (see  reference  (a)  for  example) 

t*')  X|>ti  ~  ^p+1  ^  ^  ^pff) *  Ho ' ( 4-  Ho'Gv^^X^^iy  Ho  ) 

'  (Y(f4t) "  H^ ' 

(ii^  Xp+i  ~  A  ('if4i”'tp)  •  Xp 

/s  \*T* 
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with  initial  set  up: 
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Here  only  3  by  3  matrices  need  be  inverted,  since  the  only'  inversion 
required  in  equations  Ci)-Civ)  is  of 

( 'Z  /p+n  +  (  )^  HI  )*  ) 
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Next,  at  the  end  of  m  observations  for  a  given  run,  we  obtain  X  which 
gives  the  entire  estimated  equations  of  motion:  ® 


y  ii)  =  ^x' 

'Y  ti) 

.  2(ii  =  7^y,  4 


If  we  let  fp  be  the  maximal  effective  missile  range,  then  we  consider 
the  intersection  of 

^  Sphere  of  radius  Co  centered  at  origin 

J  and 

p  . 

I  Equations  of  motion  s^,., 


(This  intersection  may  be  vacuous.) 
Thus  we  must  solve 
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the  estimated  ^notion  of  T  will  first  enter  the  maximal  missile  range, 
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The  fire  control  radar  doctrine  is  then  to  pivot  by  (random)  angles 
^  and  Tf| ,  based  on  the  data,  about  its  source  0  (see  figure  1). 

Next,  for  the  same  run,  estimated  entrance  and  exit  times  and  positions 
of  T  are  computed 

A  * 

^(m)  estimated  time 

^(m)  estimated  time 

L  1,^.2  ^0  are  predetermined  constants. 

Then  the  fire  control  radar  is  turned  on  at  time  and  is 

turned  off  at  min  ^^m)  ^^2^ 

Next,  we  consider  the  0-1  random  variable  Z  with  parameter  CO 
where 

Z  =  1  iff  the  (deterministic)  target  T  in 
following  its  actual  path  is  in 
lethal  region  L  between  on-times 
and  off- times 

=  0  otherwise 
Pr(Z=l)  =  CJ 

Thus,  for  a  given  run  we  assign  an  outcome  value  of  Z  accordingly  and^ 
best  estimate  W  by  the  percentage  of  successful  tcwls  J  that  is,  CO 
number  of  I's  over  a  period  of  runs. 

Finally,  we  let,  after  a  period  of  runs  (for  the  same  generated^ 
actual  path  of  T  and  fixed  set  of  covariance  matrices  of  errors) 
probability  of  detection  of  T  by  radar  s:  €> 


SAED  TM  70-014-3 


4S 


'  I  i  X(,^,  ^  XlM) 


•$  Irtrtpr  Mitjlf 

^0  'S  OiAcf  MhsJilc 
R«n*)e 

^•etl^il  Region  A^mas'I 

T  i»  L 

"  Es'til''>tf^  T* 

5*se<<  of\  hU  vf  to 

Tihdf  tvrt 


-4— \ 

I  6»i/rt.lM5  hr  Mjtr  ^ntti 
E*!.  l^(fosf)x+{sinviy-{t4)tri)2so 
/  ^1- I :  Cf«vix  -usinv;  y-cot<4+w).z*o 

r  Ei.ofll.'  X  -(fotvy^y  -0 

E^.c<J2:  X  =0 


FIGURE  1.  LETHAL  DETECTION  REGION  FOR  FIRE  CONTROL  RADAR 
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APPENDIX  A 

OBTAINING  OF  OBSERVATIONAL  DATA 


/.  °2'^"  '?  ‘“fl®"'"'  ‘h®  '"odel.  the  critical  parameters  r*.  e* 

trolcal  view°o^’a°  °*’'“ncd.  Figure  A-1  Illustrates ‘the^’ 

JSsent.  ■  *  operator  (azimuth-range)  when  jamming  Is 
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(See  page  4  for  other  relations  between  x  and  y  coordinates.)  ^ 

We  need  and  as  alternate  values,  direct  range  and 

See  figure  A-3. 

The  extreme  elevation  heights  of  the  target  at  time  ti  are 
O  ■<  ,  a  priori  known  beforehand.  We  let 

.'1  *  ■ 
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and  $.  #  t  (t^T)}) 
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4, 


As  before,  we  can  obtain  an  estimate  of  0^  by 

^  Arcsin  (a{‘'V 


Similarly, 
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A!''  -  A?' 


We  note  that,  in  general,  although  when  jamming  is  present,  no 
direct  visual  observation  of  the  target  can  be  made,  the  centered 
coordinates  of  the  jamming  pattern;  that  is,  the  computed 

S;  .  (0;  ,  «;  )  .  V;  ,  r,  ,  A;,  x,',  Y,,2; 

are  taken  as,  in  a  sense,  best  uniform  estimates  of  the  target’s 
coordinates  at  time  ti.  Hopefully,  biasedness  caused  by  this  procedure 
will  be  minimal. 

A  human  factors  study  is  needed  to  evaluate  for  each  i,  the  upper 
and  lower  bounds  of  the  various  parameters  - 

er.§;“-  ( 


and  relate  these  meaningfully  to  confidence  levels  for  the  corresponding 
true  parameter^  unbiasedly. 
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FIGURE  A-2.  GEOMETRY  OF  TWO  INDEPENDENT  AZIMUTH-RANGE  VIEWS 
BY  RADAR  OPERATORS  WIEN  JAMMING  IS  PRESENT- 
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FIGURE  A-3.  ELEVATION  PARAMETERS  FOR  ACQUISITION  RADAR 
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